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O Abstract 



Over the past years, the freely available Monte Carlo-code REAS which simulates radio emission from air showers based on the 
CN| geosynchrotron model, was used regularly for comparisons with data. However, it emerged that in the previous version of the code, 
emission due to the variation of the number of charged particles within an air shower was not taken into account. In the following 
QJ article, we show the implementation of these emission contributions in REAS3 by the inclusion of "end-point contributions" and 
C/D discuss the changes on the predictions of REAS obtained by this revision. The basis for describing radiation processes is an universal 
f"**) description which is gained by the use of the end-point formulation. Hence, not only pure geomagnetic radiation is simulated with 
t-H REAS3 but also radiation due to the variation of the net charge excess in the air shower, independent of the Earth's magnetic field. 
Furthermore, we present a comparison of lateral distributions of LOPES data with REAS3-simulated distributions. The comparison 
shows a good argeement between both, data and REAS3 simulations. 
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1. Introduction 

o 

£j In recent years, radio detection of cosmic ray air showers has 
55 been developed further. With radio detector arrays like LOPES 
l_ ifl], H and CODALEMA fl, fl, correlations of the radio 

signal with air shower parameters are studied and the domi- 
T 7^ nance of the geomagnetic emission contribution was verified. 

To study the physics of cosmic rays using radio signals, de- 
0^ tailed theoretical simulations are needed. Many approaches for 
0N the modelling of radio emission exist, but presently, there are 
1 J two major approaches, both of which are based on geomag- 
0^ netic effects |5J|. On the one hand the geosynchrotron model as 
<0 implemented in REAS developed by Huege et al. |6],|7],|8], 
and on the other hand the macroscopic geomagnetic radia- 
T! tion model (MGMR) of Scholten, Werner and Rusydi ifloll.lfTll. 
J> So far, the two models made conflicting predictions for the ra- 
k> dio emission of cosmic ray air showers. Essentially, this could 
5_j be seen in the different pulse shapes (unipolar for REAS2 and 
C3 bipolar for MGMR) and the differences in the frequency spectra 

for low frequencies (dropping to zero for MGMR and levelling 

off for REAS2). The details of these differences and a com- 
parison of both models are discussed in 111 211 . It arised that in 

REAS2, radiation due to the variation of the number of charged 

particles in EAS was not considered as it is the case for nearly 

all time-domain approaches as well. The reason for this missing 

contribution was a flaw in the implementation of the radiation 

process of the geosynchrotron radiation. In REAS3 [13], this 

flaw was solved. In the following sections, the details of the 

implementation as well as the results are illustrated. 
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2. General functionality of REAS 

To understand what was missing in the techincal implemen- 
tation of radio emission in REAS2 it is helpful to know the gen- 
eral structure of this Monte Carlo-code. First, the air shower is 
simulated with CORSIKA Q 

saving all important informa- 
tion, e.g. the distribution of energy and momentum of the parti- 
cles, in histograms. On the basis of these histograms, in REAS, 
shower particles are generated according to the desired distri- 
butions derived with CORSIKA. In the simulation code, each 
particle is followed analytically on its track through the Earth's 
magnetic field. Note that the real particle trajectories are de- 
scribed by several unrelated short tracks. Finally, the radiation 
given from all shower particles is superposed for each single 
observer position. In REAS2, only radiation processes along 
the trajectories were treated, but not at the end or the begin- 
ning of the tracks. This can be compared with a situation that 
the particles arrive with velocity v m c given by CORSIKA, 
enter the Earth's magnetic field where they are deflected on a 
short curved track and finally fly out of the influence of the ge- 
omgagnetic field with velocity v ~ c (cf. left sketch of Fig. 
[TJ. To revise the flaw in the derivation of the radiation process 
of the geosynchrotron radiation, emission contributions at the 
beginning and the end of the tracks have to be taken into ac- 
count, i.e. radiation due to the acceleration of the particle at the 
starting point of the trajectory and vice-versa deceleration of 
the particle at the stopping point is considered. To implement 
this radiation, the best description of the radiation processes is 
an end-point formulation 11511 . In this case, the tracks are de- 
scribed by straight track segments joined by "kinks" (cf. right 
sketch of Fig. [TJ. If at a given atmospheric depth more particle 
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Figure 1: Sketch of the trajectories how they are implemented in REAS. Left: 
REAS2. Right: REAS 3. 



trajectories start than end, e.g. the number of particle declines, 
this results in a net contribution. 



3. Incorporation of end-point contributions 

Adding the discrete end-point contributions to the continu- 
ous contributions along the tracks may produce problems, e.g., 
there is a risk of double-counting. To get a consistent descrip- 
tion of all radiation processes in the simulation, it is convenient 
to use the end-point formulation throughout. Radiation occurs 
if the velocity of the particle changes, i.e., in a kink of the track. 
Because the change of the velocity can be considered instanta- 
neous with respect to the times of interest (6t <k — - — and 

r "observed 

Vowed < 100 — 1000 MHz) only the time-averaged process is of 
interest. Hence, the time-integrated field strength of the radia- 
tion formula can be calculated. Equation[TJshows the result for 
the radiation in one kink of the track, 
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a simple shower geometry is chosen where the geomganetic 
angle is 90°, i.e., a vertical shower with a primary energy of 
10 17 eV and a horizontal magnetic field of 0.23 Gauss was se- 
lected. Since one typical shower out of many CORSIKA sim- 
ulated air showers was chosen, shower-to-shower fluctuations 
do not influence this comparison. For REAS2 and REAS3 the 
same CORSIKA shower was taken as the basis. In Figf2]the 
raw pulses of REAS2 and REAS3 for an observer 100 m north 
of the shower core are shown as well as the frequency spectra 
for observers 100 m north and east of the shower core. It is obvi- 
ous that the pulse shape changed from unipolar to bipolar. This 
change agrees with the theoretical expectation since the source 
of the radio emission exists only over a finite time in a finite re- 
gion of space (cf. HH). ln me frequency spectra (right plot of 
Fig. |2]i this behaviour can be seen as well because the spectral 
field strengths drop to zero for frequency zero. In addition, the 
spectral field strength for observers at different azimuthal posi- 
tions differs less in REAS3 than in REAS2, which indicates an 
increased azimuthal symmetry of REAS3 compared to REAS2. 
For the spectra of Fig. [2] one observer in the north and in the 
east of the shower core is selected. In general, the spectra got 
flatter for REAS3. The increased azimuthal symmetry is again 
visible in the contour plots of Fig. [3]for the total field strength 
which gives an overall impression of the changes from REAS2 
to REAS3. In the contour plot of REAS3, an east-west asym- 
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Figure 3 : Contour plots of the 60 MHz absolute field strength for the radio emis- 
sion from a vertical air shower. Left: REAS2. Right: REAS3. The white region 
in the center has not been simulated. Contour levels are 0.1 /iVm~'MHz~' 



where e indicates the particle charge, p = $(t)/c is given by the 
particle velocity, R(t) = \R(t)\ describes the vector between par- 
ticle and observer position and n{t) - R{t)/R(t) is the line-of- 
sight direction between particle and observer. The index "ret" 
means that the equation needs to be evaluated in retarded time. 
fi\ corresponds to the velocity before and to the velocity af- 
ter the kink. In this completly universal (cf. 1151 ') and discrete 
calculation, radiation at the end or the beginning of the track 
corresponds to kinks where one velocity is equivalent to zero, 
i.e., one term of the integrated sum vanishes. 

4. Results 

4.1. Comparison between REAS2 and REAS3 

In this section, a short overview about the major changes 
from REAS2 to REAS3 is given. In II 1311 more details and dis- 
cussion can be found. For the comparison, several simulations 
were done with a set of prototype showers. For this article, 



metry is distinguishable, i.e., the signal in the east is higher than 
the signal in the west. This is explainable with the fact that in 
an air shower a charge excess of electrons occurs. The variation 
of the net charge excess leads to a further emission contribution 
which is discussed in the following section. 

4.2. Discussion of charge excess 

Due to the time variation of the net charge excess in air show- 
ers, radio emission occurs even in the absence of any magnetic 
field. The same effect appears if the incoming direction of the 
air shower is parallel to the Earth's magnetic field. In REAS2, 
there was no radio emission for both cases because radio emis- 
sion arised only due to the deflection in the magnetic field. For 
a pure emission due to variation of the net charge excess a radi- 
ally polarised component is expected. To verify this, the radio 
emission of a vertical air shower with primary energy of 10 17 eV 
was simulated while the strength of the magnetic field was set 
to Gauss. The contour plots of the 60 MHz field strength for 
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Figure 2: Direct comparison of REAS2 (dashed blue) and REAS3 (solid red) for a vertical air shower and observer distance of 100 m. Left: Raw pulse for an 
observer 100m north of the shower core. Right: Frequency spectra for observers east (thick lines) and north (thin lines) of the shower core. 



the north and the west polarization components are shown in 
Fig. [4] The emission pattern is, as expected, radially polarised 
but the relative field strength of the charge excess emission at 
60 MHz is small in the distance range up to 200 m compared to 
the relative field strength of total radio emission given in Fig[3] 
(the contour levels for the simulation in this section are a factor 
of ~3.3 smaller than for the simulation with a realistic magnetic 
field). Due to the charge excess in air showers and its varia- 
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Figure 4: Contour plots of the 60 MHz field strength for the charge excess 
component from a vertical air shower. Left: north polarization component. 
Right: west polarization component. The white region in the center has not 
been simulated. Contour levels are 0.03 /iVnr'MHz" 1 

tion it is evident that in the footprint of the radio signal a small 
azimuthal asymmetry has to remain. Since the charge excess 
leads to a non-geomagnetic contribution in the radio emission, 
it is obvious that the radio emission of extensive air showers is 
not purley v x B dependent. 

4.3. Comparison with LOPES data 

Finally, to test the predictive capability of the REAS3 simula- 
tion on real data, it is suitable to compare the output of REAS3 
and REAS2 with measured data. In this case, the lateral dis- 
tributions of LOPES data have been used (published in lfl7lo . 
First, a set of 200 showers per each event was simulated with 
CONEX. Then, a "typical" shower was chosen out of this set 
to avoid shower-to-shower fluctuations. For the selection, the 
averaged depth of the shower maximum X max of all 200 show- 
ers was calculated and then a shower was chosen with an X max 



similar to the averaged one. For REAS2 and REAS3 the same 
CORSIKA shower was used to simulate the radio emission. In 
Fig f5] two example events are shown for showers with proton 
as a primary particle. The LOPES event shown in the left plot 
has a primary energy of 2.9 • 10 17 eV and a zenith angle of 31°. 
The REAS3 simulation reproduces the measured data very well. 
The second event has a primary energy of 2.9 • 10 18 eV while the 
zenith angle is 58°. For the first time, it is possible to reproduce 
also some flat lateral distributions which account for ~ 10% 
of all measured events. In general, the lateral distribution with 
REAS3 became flatter than with REAS2 and agrees much bet- 
ter with the data. There are still some events which are not re- 
produced correctly by REAS3 but for all events the simulations 
are much closer to the measured data than REAS2 simulations. 
For future investigations, the selection of a typical shower will 
be changed such that the ratio of the electron and muon num- 
ber is comparable to the number given by KASCADE instead 
of the selection via the X max . This could increase the agreement 
between simulations and measurements as well as the noise cor- 
rection lfl8ll on the data. For further comparisons between sim- 
ulations and data a detector simulation, which is not included 
so far, should be implemented to take into account the com- 
plete influence of the detector characteristics on the measured 
signal. Since REAS3 (the same is true for REAS2) has no free 
parameters, the simulations are only dependent on the input pa- 
rameters given by the shower characteristics such as primary 
mass, energy and incoming direction. 

5. Conclusions 

The incorporation of end-point contributions as described 
in the article, led to a revision of the implementation of the 
geosynchrotron model in the Monte Carlo-code REAS3. Due 
to the universal end-point formalism not only emission contri- 
butions due to the variation of the number of charged particles 
in an air shower are taken into account but also contributions 
due to the deflection of the particles in the geomagnetic field as 
well as a radiation component arising due to the variation of the 
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LOPES event 1 161 186603 filtered lateral 
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Figure 5: Lateral 
showers. 



distributions of two different LOPES events (blue *) compared with REAS2 (pink +) and REAS3 (red x) simulations of proton induced air 



net charge excess in an extensive air shower. With this revision, 
the pulse shape became bipolar and the frequency spectra are 
dropping to zero for zero frequency. Remaining asymmetries 
in the nearly symmetric azimuthal emission pattern are explain- 
able by the time-varying charge excess in air showers. This re- 
sults in a radio emission which is not purely of geomagnetic ori- 
gin. REAS3 shows a good agreement between simulation and 
LOPES data comparing the lateral slopes, which is particularly 
striking since RE AS 3 has no free parameters at all. RE AS 3 
is the first self-consistent time-domain model which takes the 
full complexity of air shower physics as provided by CORSIKA 
into account. The code will be freely available for users world- 
wide. 
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